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Abstract 



A search for charginos with masses close to the mass of the lightest neutralino is 
reported, based on the data collected with the DELPHI detector at LEP from 
1995 to 1997 at centre-of-mass energies between 130 and 183 GeV. The signature 
of a photon at high transverse momentum radiated from the initial state reduces 
the two-photon background to acceptable rates, thus making the mass differ- 
ences between a few hundred MeV/c 2 and 3 GeV/c 2 detectable. In very nearly 
degenerate scenarios, the lifetime of the chargino can be large enough to produce 
either visible secondary vertices or decays outside the detector; therefore, quasi- 
stable heavy charged particles and displaced decay vertices were also searched 
for. No excess of events with respect to the Standard Model expectations was 
observed, and limits in the plane of chargino-neutralino mass difference versus 
chargino mass are given. 
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1 Introduction 

Supersymmetry (SUSY) [|IJ is an appealing theory which answers some well-known 
outstanding questions of the Standard Model (SM), at the expense of introducing super- 
symmetric partners of the known particles (sparticles). The most well-studied example 
of supersymmetry is the minimal supersymmetric extension of the SM (the MSSM). If 
the .R-parity quantum number is conserved, as often assumed, there must exist a lightest 
supersymmetric particle (LSP) which is stable and remains after any SUSY decay chain. 
Such an LSP is expected to be neutral and weakly interacting ||. The usual way of 
searching for SUSY particles heavier than the LSP in e + e~ interactions is therefore to 
look for visible particles accompanied by the missing energy carried away by two or more 
LSP's. This works whenever the mass difference between the produced sparticle and the 
LSP is large enough to leave some sizeable amount of energy for the visible final state 
particles. Typically, the searches carried out so far at LEP in the different SUSY chan- 
nels go down to mass differences of a few GeV/c 2 . Previous DELPHI searches at LEP2 
have resulted in limits on the production of charginos, neutralinos, sleptons and b and t 
squarks in the MSSM, valid when the mass difference between the sparticles searched for 
and the LSP (usually the lightest neutralino) is above 3 to 5 GeV/c 2 |||]]. For smaller 
mass differences, the only limits available so far are those derived from the precise mea- 
surement of the Z width at LEP1; in particular charginos lighter than about Mz/2 are 
excluded, irrespective of their field content |||J. 

The search for charginos and neutralinos is essential to constrain SUSY. It is therefore 
of paramount importance to make sure that they have not been missed because of possible 
small mass differences amongst them. Such small mass differences are rather unlikely 
in the MSSM if the masses of the gauginos are assumed to be all the same at some 
grand unified scale, as expected in supergravity (SUGRA) models. In such models the 
lightest chargino and the two lightest neutralinos can have nearly the same mass only if 
those gaugino masses are unnaturally large (above 1 TeV/c 2 ). However, since no direct 
experimental support for those models has been found so far, it is reasonable to also 
consider models without the SUGRA assumptions. In particular, there are interesting 
theoretical string-motivated scenarios which explicitly prefer the non-unification of the 
gaugino masses at the GUT scale |7] and in which is quite likely that the lightest chargino 
and the lightest neutralino have nearly equal masses ||. 

The region of low mass difference is experimentally challenging. If the chargino- 
neutralino mass difference AM ± = M-± - o is below the mass of the pion, the lifetime 
of the chargino can be so long that it passes through the entire detector before decaying. 
In DELPHI this can be covered by a search for heavy charged particles identified with 
the Cherenkov detectors and/or because of their anomalously high ionization in the gas 
chambers. Mass differences of a few hundred MeV/c 2 may be observed by looking for re- 
constructed secondary vertices from a chargino decay inside the detector but significantly 
displaced from the main interaction point. 

When AM increases, the chargino decay length becomes so short that the decay 
vertex can hardly be distinguished from the production one. However, as long as AM^ 
remains below a few GeV/c 2 , the visible particles carry only a small fraction of the parent 
energy. A minimum visible transverse momentum is usually required to reject two-photon 
interactions. This may result in an almost complete loss of efficiency for chargino decays. 
On the other hand, some transverse momentum requirement is necessary because the two- 
photon cross-section is orders of magnitude higher than any signal searched for at LEP2. 
Here a suggestion by Chen, Drees and Gunion H is applied to search for these charginos 
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at low AM 1 * 1 . If one considers the events accompanied by a hard photon from Initial 
State Radiation (ISR), then the two-photon background can be kept small by choosing 
the photon transverse energy to be greater than 



1 + sin 9 min 

where 9 m i n is the lowest polar angle accessible in the detector. If an ISR photon with 
a transverse energy above (K^) min is radiated from a two-photon event then, typically, 
one of the final state electrons, which usually escape undetected in the beam pipe, should 
be deflected at an angle larger than 9 m [ n , thus allowing the identification of the event 
as background. Such a selection gives a low efficiency since only a small fraction of the 
SUSY events have an ISR photon with ET > (E^) mm . On the other hand, the presence 
of the high energy photon in the detector substantially increases the otherwise low trigger 
efficiency for these decays with only a few soft visible particles. 

This paper first explores the feasibility of a search at LEP2 for charginos or second 
lightest neutralinos nearly mass-degenerate with the lightest neutralino. In the case of the 
second lightest neutralino (for which the relevant mass difference is AM = M^o — M^o) it 
will be shown that a search at LEP2 is either impossible or quite difficult, at least for most 
of the values of AM of interest here. Instead, a search for mass-degenerate charginos was 
found to be feasible, and was realized using the data collected by the DELPHI experiment. 
For sensitivity in the case of long chargino lifetimes two alternative searches were used, 
as described in section |j: one is the search for heavy stable charged particles described 
in ||10| , and the other is a modified version of a search for secondary vertices |11| with 



reconstructed incoming and outgoing tracks (kinks). The search which exploits the ISR 
signature to cover the mass differences between 0.3 and 3 GeV/c 2 was specifically designed 
for this work and is discussed in section |^. 



2 Data samples and event generators 



A detailed description of the working DELPHI detector can be found in JT2[. The 
trajectories of the charged particles are reconstructed in the 1.2 T magnetic field by a 
system of cylindrical tracking chambers. The most relevant for the analyses reported 
here are the Silicon Tracker, the Inner Detector (ID), and the Time Projection Chamber 
(TPC). The Silicon Tracker is composed of the Vertex Detector (VD) in the barrel and 
the ministrips and pixels of the Very Forward Tracker (VFT) at low polar angles (#); it 
covers the range between 10° and 170° in 9 and radii down to 6.3 cm from the beam. 
The ID covers polar angles down to 15° (165°). The TPC tracks particles between the 
radii of 29 and 122 cm, with at least three pad rows crossed if 9 is between 20° and 160°. 
The electromagnetic calorimeters are the High density Projection Chamber (HPC) in the 
barrel (40° < 9 < 140°) , the Forward ElectroMagnetic Calorimeter (FEMC) in the forward 
regions (11° < 9 < 36° and 144° < 9 < 169°) and the Small angle Tile Calorimeter (STIC) 
in the very forward part (down to 1.66° from the beam axis). In front of the STIC, which 
is also the luminometer of DELPHI, have been placed two planes of scintillators (Veto 
Counters |TB[) used to detect charged particles which enter the calorimeter. Excellent 
particle identification is provided by the Ring Imaging CHerenkov (RICH) detectors, 
equipped with two different radiators (liquid and gas) with different refractive index, and 
thus different momentum thresholds. 

During the high energy runs of LEP in 1995-97, DELPHI collected data at centre-of- 
mass energies of approximately 130, 136, 161, 172 and 183 GeV. Only the runs in which 
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the relevant subdetectors worked correctly were taken into account for each analysis. The 
luminosities used at the different energies in the analyses considered here are approxi- 
mately 6, 6, 10, 10 and 54 pb _1 respectively. The data at 130 and 136 GeV where not 
searched for kinks; in the search for soft particles accompanied by ISR only 50 pb _1 of 
the run at 183 GeV could be used, mainly because of some temporary degradation of the 
quality of the data collected by the HPC and FEMC calorimeters, which are fundamental 
for that analysis. 

To evaluate the signal efficiency and the background contamination, events were gen- 
erated using several different programs, all relying on JETSET 7.4 |14| for quark frag- 
mentation. All the events generated were passed through a complete simulation of the 



DELPHI detector fllS] , and then processed in exactly the same way as the real data. 

The program SUSYGEN ||16|| , which includes initial and final state photon radiation, 
was used to simulate all signal events. The implementation of the decay widths (i.e. 
branching ratios and lifetimes) in SUSYGEN at low AM ± has been modified in order to 



reproduce the results of the analytical calculations reported in Jy and JL7]] . 

For the Standard Model background, several samples of the different final states were 
generated with statistics which were typically well above those expected (although some 
of the two-photon samples, especially at the lowest centre-of-mass energies studied here, 
were originally simulated with statistics only slightly higher than the one expected). An- 
nihiliation of e + e~ into a virtual Z/7, including ISR, were generated with PYTHIA fT4|| . 
This generator was also used for four-fermion processes. For the two-photon collisions, 
the generator of Berends, Daverveldt and Kleiss (BDK) [|18j was used for the leptonic final 
states. In this generator the (e + e~)/i + /i~ and (e + e - )r + r~ final states include the simula- 
tion of the transverse momentum of the ISR photon, while in (e + e~)e + e~ only collinear 
ISR is implemented. In the two-photon interactions leading to hadronic final states, the 



QCD and VDM parts were simulated with the TWOGAM [T9j generator in which ISR 



is generated without transverse momentum. BDK, with visible ISR, was used for the 
QPM part. Although not used for the computation of the background, QPM events were 
also generated using the TWOGAM generator, to study the differences between the two 
programs, in particular effects related to the transverse momentum of the ISR photons. 
The absence of these effects in some of the available two-photon samples, implies that the 
simulation underestimates the background from ISR in two-photon events. 



3 Chargino and neutralino production and decay at 
low AM 

The higgsino and gaugino sector of the MSSM can be described in terms of four pa- 
rameters: the ratio tan f3 of the two Higgs vacuum expectation values, the Higgs mixing 
parameter /i, the SU(2) gaugino mass M2 and the U(l) gaugino mass Mi. In the models 
with gaugino mass unification at the GUT scale, there is a relation between Mi and M2 

5 

Mi = - tan 2 9 W ■ M 2 ~ 0.5 • M 2 . (2) 
3 

However, it has been already anticipated that this unification is not strictly necessary in 
the theory and there are several models without it, in particular the SUSY-string scenario 
proposed in ||. The definition 



Mi = R f ■ - tan 2 6 W ■ M 2 
3 



(3) 
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will be used, so that any value of Rf different from 1 indicates how much the model 
deviates from the gaugino mass unification hypothesis. 

In the MSSM there are two charginos (xf and xf). These mass eigenstates are linear 
combinations of the two interaction eigenstates, the wino and the charged higgsino. There 
are also four neutralinos, linear combinations of the neutral interaction eigenstates. In 
the following it will be assumed that the lightest neutralino is the LSP. 

The lightest chargino gets almost the same mass as such an LSP in two cases ||: 

1. Low large M 12 scenario: Xi an d Xi are both higgsino-like and nearly degenerate, 
with masses ~ 

2. High |/i|, low M2 scenario: the Xi an d the xf are both gaugino-like and nearly 
degenerate with masses ~ M%. In this scenario, in order to have AM around 1 
GeV/c 2 or smaller, Rf in eq. (§) must be larger than or equal to 2. 

In the first scenario, the second lightest neutralino is also almost mass-degenerate with the 
lightest neutralino, with a mass splitting which is sligtly larger than that of the lightest 
chargino. The same is true also for the second scenario, but only for Rf ~ 2. 



3.1 Cross-sections 



The predicted e + e~ — > XiXi an d e + e~ — > X1X2 cross-sections |L6| at the centre-of- 
mass energy of 183 GeV are shown in figures |l| and ||, respectively, as functions of the xf 
and X2 masses. That value of the centre-of-mass energy is taken as example, since the 
behaviour is similar for all energies studied. In both figures, the upper plot refers to the 
higgsino cross-section and the lower plot to the gaugino one. 

For the neutralinos, the e + e~ — > X1X1 an d e + e~ — > X2X2 cross-sections are much smaller 
than the e + e~ — > X1X2 one 5 shown in figure |2|. This last process is the only one considered 
in the following for the production of neutralinos at LEP2. 

The widths of the bands arise from a variation of tan/3 from 1 to 50, M2 and |//| 
below 100 TeV/c 2 , Rf between 1 and 10. The higgsino cross-sections are quite stable 
when varying Rf, and also the values assumed by the charged gaugino cross- sections do 
not deviate by more than 5% if Rf moves in that range. In the chargino production, 
the exchange of sfermions in the t-channel interferes destructively with the s-channel 
production, lowering the gaugino cross-section at small values of M$. On the contrary, in 
the neutralino production the interference is constructive and the gaugino cross-sections 
are enhanced at low Mg. In all the figures the mass of the relevant scalar lepton has been 
varied from to 1 TeV/c 2 . 

The chargino cross-section, in the approximation of large scalar masses, is roughly three 
times larger in the gaugino-like scenario than in the higgsino-like one. On the contrary, 
the e + e~ — > X1X2 cross-section in the gaugino-like scenario is typically several orders 
of magnitude smaller than in the higgsino one, certainly out of reach of the luminosity 
planned at LEP2. 



3.2 Branching ratios and lifetimes 

The partial decay widths of the chargino in SUSYGEN have been modified in order to 
account for the decays into a neutralino and one, two or three pions, according to analytical 
calculations ||[17||. These calculations use the form factors of the low mass hadronic 
resonances to determine the hadronic width for < AM^ < 2 GeV/c 2 . This treatment 
has been verified || to describe correctly the r decays, i.e. decays with AM = m T . 
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Figure 1: Predicted e + e~ — > XiXi cross-sections at the centre-of-mass energy of 183 
GeV as a function of the mass of the chargino. The upper plot refers to the -C M 2 
higgsino-like scenario; the lower plot to the M 2 <C \fi\ gaugino-like scenario. The widths 
of the bands allow for a variation of Mi, M 2 and fi so that < AM < 3 GeV/c 2 ; 
1 < tan/5 < 50; M 2 < 2M 1 < 10M 2 ; Afc > M~± (the upper part of the gaugino band 

displays separately the points corresponding to My > 500 GeV/c 2 ). 



Figure |]a and b show the leptonic and hadronic branching ratios computed for a 
50 GeV/c 2 charged higgsino with a mass between 100 MeV/c 2 and 5 GeV/c 2 above that 
of the lightest neutralino. In the plot of the hadronic modes, the decays into a Xi an d one, 
two or three n, which contribute to the total xt ~^ Xi<¥l width, are shown separately For 
a gaugino, the leptonic width is enhanced for low M„ because of the sneutrino mediated 
decays, and this was taken into account in the analysis. 

Figure |3|c shows the lifetime of the same chargino as a function of AM 1 . The figure 
clearly shows the step caused by the onset of the dominating xt X? 7r+ t wo body decay. 
The chargino lifetime and branching ratios depend strongly on AM 1 and relatively little 
on the SUSY scenario, as long as scalar exchange can be neglected. For low Mp, charged 
gauginos get a shorter lifetime, and this was also taken into account in the analysis. 

Figure f| shows the exclusive branching ratios of the second lightest neutralino as a 
function of AM , in the two scenarios which allow almost degenerate states. Unlike the 
chargino case, there is now a strong model dependence of the decay widths. As already 
mentioned, the only neutralino scenario with a sufficient cross-section to be searched for is 
the higgsino one. In this scenario the rate of the radiative decay x® ~ * X?7 increases in the 
low AM region of interest here. Thus a large fraction of the x°> decays yield either a low 
energy photon or a pair of neutrinos, in addition to the LSP. Since it is difficult to identify 
photons below about 1 GeV in the detector and distinguish them from background, both 
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Figure 2: Predicted e + e~ — > X1X2 cross-sections at the centre-of-mass energy of 183 
GeV and as a function of the X2 mass. The upper plot refers to \fi\ <C M 2 higgsino-like 
scenario; the lower plot refers to the M 2 <C gaugino like scenario. The width of the 
bands allows for a variation of M2, Mi and [i so that for any given mass of the second 
neutralino < AM < 3 GeV/c 2 ; 1 < tan/5 < 50; 1 < R f < 10 (R f ~ 2 for the gaugino); 



decay modes can be considered invisible for practical purposes. For these reasons the 
present work has been limited to charginos, leaving neutralinos aside. 



4 Search for long-lived charginos 

Two methods were used to look for charginos with a visible decay length: the search 
for heavy stable charged particles and the search for decay vertices inside the detector. 
Both searches are described in this section. 



4.1 Heavy stable charged particles 

Heavy stable or almost stable charged particles are identified through their anomalously 
high specific ionization in the TPC or by the absence of Cherenkov light produced in the 



two radiators of the barrel RICH. The leptonic selection described in |10[ is the one used 
for the present analysis. The efficiency of this selection for pair-produced heavy charged 
particles, traversing the full depth of the detector, is given in [|10| . Figure || shows the 



efficiency for selecting a single heavy charged particle as a function of its mass and of the 
LEP centre-of-mass energy. The trigger efficiency for high momentum charged particles 
crossing the full depth of the TPC and the RICH is practically unity. 
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Figure 3: Predicted (a) leptonic and (b) hadronic branching ratios and (c) lifetime of a 
50 GeV/c 2 chargino (higgsino) decaying into a Xi P ms standard particles, as a function 
of the mass difference between the xf an d the Xv Similar decay modes are predicted 
for the gaugino scenario in the approximation of large M„, which makes the contribution 
of the v mediated decays negligible. For smaller M$, the leptonic decays of the charged 
gaugino are enhanced and the lifetime gets shorter. 




Figure 4: Predicted branching ratios of the decay modes of the second neutralino, with 
M^o ~ 50 GeV/c 2 , as a function of the mass difference between the \2 an d the Xi- The 
left and right plot refer to the higgsino-like and gaugino-like scenario respectively. 



Lower limits on the mass of the chargino have been already published [TO] under the 
hypothesis that it decays predominantly outside the DELPHI detector (in case of an heavy 
v this corresponds to AM < 100 MeV/c 2 ). For shorter lifetimes, the detection efficiency 
can be obtained by convoluting the efficiency for stable particles with the probability that 
the chargino passes through the barrel RICH before decaying. 



4.2 Decay vertices inside the detector (kinks) 

If the heavy charged particle decays inside the central tracking devices of DELPHI (at 
a radius between 10 cm and 1 m) then both the incoming and the outgoing track can be 
reconstructed, and the angle between the tracks can be calculated. This method was used 
in a DELPHI search for scalar tau leptons decaying into a light gravitino and an ordinary 
tau lepton For small AM ± , however, the visible momentum of the decay products 
is quite small (typically less than 1 GeV/c), and the identification of the secondary track 



and, therefore, of the kink becomes more difficult. The selection criteria adopted in [[□]] 
have therefore been modified for the present search, exploiting the typical topology of 
these events: two particles emitted in opposite hemispheres decaying into one low-energy 
charged particle each. 

A set of rather loose general requirements was imposed on the events in order to 
suppress the low energy background (beam-gas, beam- wall, etc), two-photon, e + e~ and 
hadronic events: 

• there must be at least one charged particle and not more than five; 

• the visible energy must be above 10 GeV; 

• the total energy in electromagnetic showers was required to be below 60 GeV; 

• the transverse momentum with respect to the beam axis had to be greater than 
5 GeV/c; 
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Figure 5: Efficiency for detecting a single heavy stable chargino in DELPHI (produced in 
pair in the e + e~ rest frame), as a function of its mass and for the different centre-of-mass 
energies used in the analysis. 



• the energy measured in the STIC must not exceed 10 GeV. 

To compute the above quantities the reconstructed charged particles were required to have 
momenta above 100 MeV/c and impact parameters below 4 cm in the transverse plane 
and below 10 cm in the longitudinal direction (however, no requirement on the impact 
parameters were imposed to the reconstructed tracks in the following steps). Clusters 
in the calorimeters were interpreted as neutral particles if they were not associated to 
charged particles and if their energy exceeded 100 MeV. 

All the charged particle tracks were grouped in clusters according to their measured 
point closest to the interaction vertex (starting point). The clustering procedure is de- 
scribed in [[Ll|]. Each cluster contains all tracks whose starting points differ by less than 
2 cm. The starting point of a cluster is defined as the average of the starting points of its 
tracks. This procedure allows for clusters with a single track. 

A cluster with only one track was considered a Xi candidate if: 

• the distance of the starting point from the beam spot, in the plane transverse to the 

beam axis (xy plane), R S p , was smaller than 10 cm; 

• its momentum was greater than 20 GeV/c; 

• the polar angle of the momentum had to satisfy | cos6*| < 0.8; 

• the impact parameters of the track along the beam axis and in the plane perpendic- 
ular to it were less than 10 and 4 cm, respectively 



For each single track cluster fulfilling the above conditions, a search was made for a 

second cluster with starting point in the transverse plane beyond R S p , and an angular 
separation between the directions defined by the beam spot and the starting points of 
the clusters smaller than 90° in the xy plane. This secondary cluster was assumed to be 
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formed by the decay products of the xt ■ Therefore, the xt candidate and the secondary 
cluster had to define a vertex. If the secondary cluster included more than one track, only 
the track with the highest momentum was used to search for the decay vertex or kink 
(crossing point with the xt track). 

The crossing point was defined as the midpoint of the line segment connecting the 
points of closest approach in the xy plane between the two (possibly extrapolated) tracks: 
the candidate xt track and the selected track from the secondary cluster. The following 
conditions were required to define a good crossing point: 

• the minimum distance between the tracks had to be smaller than 1 mm in the xy 
plane, 

• the crossing point, the end point of the xt track and the starting point of its decay 
products were required to satisfy the following conditions: 

-10 cm < (R CIOSS - R^ d ) < 25 cm 
-25 cm < (R CIOSS - R^ c prod -) < 10 cm 

where R^ d , -Rcross and i?gp C - prod - are the distance from the beam spot to the end point 
of the xt track, the crossing point of the tracks and the starting point of the tracks 
supposed to come from the decay products of the chargino, in the xy plane. 

Reconstructed secondary vertices could also be the result of particles interacting in 
the detector material, or bremsstrahlung, giving a particle trajectory reconstructed in 
two separate track segments. To eliminate this kind of background, events with a good 
crossing point (kink) were subjected to additional requirements: 

• to reject hadronic interactions, any secondary vertex reconstructed in the region of 
the detector where there is material must be outside a cone of half opening angle of 
5°, with apex at the beam spot and centred around the kink direction; 

• to reject photon radiation, in the case of secondary clusters with only one track, 
no neutral particle was allowed in a 1° cone around the direction of the missing 
momentum, defined by the difference between the momentum of the xt an d that of 
the daughter; 

• to reject segmented tracks, the angle between the tracks used to define a vertex, 
calculated at the crossing point, had to be larger than 2°. 

Finally, for an event to be accepted, at least one charged particle must be found in 
each hemisphere (defined by the plane which contains the beam spot and is perpendicular 
to the line connecting the beam spot to the kink). 

This selection was applied to samples of e + e~ — > xtxl events, generated at yjs = 161, 
172 and 183 GeV and passed through the full DELPHI simulation and reconstruction 
chain. The efficiencies for the single vertex reconstruction as a function of the radial 
distance from the decay to the beam spot are plotted in figure B|a,, where the case of a 
65 GeV/c 2 chargino with AM ~ 150 MeV/c 2 has been taken as example. The selec- 
tion efficiency (e scl ) is almost independent of the decay radius for radii between 30 and 
90 cm for almost all the masses generated; it tends to decrease only when the mass of 
the chargino approaches the kinematical limit, because the momenta of the secondaries 
are not sufficiently enhanced by the reduced boost. The dependence of the selection ef- 
ficiency on AM ± is weak for the values of AM^ within the range searched for with this 
method. There is however some increase of the efficiency with increased AM , as the 
mean momentum of the decay products gets higher. 
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Figure 6: Efficiency for the kink reconstruction of a 65 GeV/c 2 chargino with AM ± ~ 150 
MeV/c 2 , as a function of the chargino decay radius, at the centre-of-mass energy of 183 
GeV: (a) selection efficiency for the single arm; (b) single particle trigger efficiency for 
the selected kinks. The lines superimposed represent the result of the fit described in the 
text. 

The trigger efficiency (e trs ) depends on the mass of the chargino and on the centre- 
of-mass energy. It has been estimated with a simplified simulation of the single track 
trigger in DELPHI, taking the overall trigger efficiency as the logical OR of the trigger 
efficiencies of every single charged track in the event [OT. This is a conservative estimate 
because the possible contribution of all higher multiplicity triggers was neglected. The 
efficiency for triggering on the hemisphere of the kink, as a function of the decay radius 
of the chargino, is shown in figure Pd, for the same sample as in figure |]a. The radial 
dependence of both selection and trigger efficiencies was fitted with polynomials, allowing 
a different fit at every mass generated. The results of these fits were used in the analysis 
and the efficiencies for masses other than those simulated were obtained by interpolation. 

4.3 Global efficiency for charginos with visible decay length 

The selections in the search for long-lived charginos yield an efficiency e(x) which is 
a function of the decay radius x of the chargino. Each chargino produced in the e + e~ 
collision can be selected by either one of the two searches: if it decays outside the RICH, 
by the search for stable charged particles; if it decays inside the ID or the TPC, by the 
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search for kinks. The global efficiency, which is a function of the decay radius, is given by 
the logical OR of the two selections. 

To extrapolate the efficiency of the event selection to the points in the space of the 
SUSY parameters not fully simulated, a semi-analytic calculation was used. The distri- 
bution of the decay length of the chargino in the reference frame of DELPHI was derived 
analytically in the scenarios studied and for any given chargino mass, AM ± and centre-of- 
mass energy. These distributions were convoluted with the efficiencies of the experimental 
search methods e(x), as determined for the fully simulated events, giving the detection 
efficiency in any point of the space of the SUSY parameters. 

As an example, figure [7] shows the combined detection and trigger efficiency at 183 GeV 
for e + e~ — > xfxi, when the charginos are pure gauginos and M v = 1 TeV/c 2 , My, = 65 
or 80 GeV/c 2 , as a function of AM ± . The efficiencies displayed include both searches for 
long-lived charginos (for heavy stable particles and for kinks). 

DELPHI 




AM (GeV/c 2 ) AM (GeV/c 2 ) 



Figure 7: Total efficiency for detecting a long-lived charged gaugino in DELPHI using the 
analysis at 183 GeV as function of AM^ (that is of its lifetime) and in the approximation 
of heavy sneutrinos. The two masses of 65 and 80 GeV/c 2 have been chosen as examples. 



4.4 Results 

The search for heavy stable particles was performed with all the data collected at 
y/s = 130/136, 161, 172 and 183 GeV. No candidates remained in the data, while 0.7±0.3 
background events were expected [ PHI . 

In the search for kinks, only data taken in 1996 and 1997 at the centre-of-mass energies 
of 161, 172 and 183 GeV were analysed. No events were selected, while the background 
expected at the three centre-of-mass energies was, respectively, 0.11 ±0.11 , 0.04^;^ and 
0.21 ±0.14 events. 

In the absence of candidates selected in any of the searches, the 95% confidence level 
(CL) upper limit is 3.0 events [ET] for the whole statistics analysed. When this limit 



is considered together with the integrated luminosities, the expected cross-sections and 
efficiencies for different SUSY parameters, the 95% CL exclusion regions in the plane 
(M-+, AM ± ) shown in fi gure [12| can be derived. As far as the SUSY scenarios are con- 



cerned, three cases were considered (see also section |||; here the results for gauginos are 
further subdivided according to the mass of the sneutrino): 
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1. low |/i|, large M 12 and any mass for the sfermions, (higgsino-like); 

2. high \n\, low M 2 and M„ > 500 GeV/c 2 (gaugino-like); 

3. high \fj,\, low M 2 and M-+ < Mj> < 500 GeV/c 2 (gaugino-like). 



i 

2 



In all cases, tan (3 was varied between 1 and 50, M 2 and between and 100 TeV/c 
and i?/ between 1 and 10 (although for Rf < 2 it is not possible to have charged gauginos 
with masses close to that of the lightest neutralino). The lowest value of the cross-section 
obtained in the scan of the SUSY parameters for each scenario, at a given mass of the 
chargino, was used for the calculation, so that the limit obtained remains valid for all 
values of the parameters. A coarse scan extending to larger values of Rf suggests that 
all the results remain valid also for Rf > 10. The limits obtained in the third case are 
much less stringent than the other two, because a light sneutrino increases the chargino 
decay width through virtual v exchange. Whenever M„ < M-+ the exchanged sneutrino 
becomes real and the lifetime drops: this region cannot be excluded by the search for 
long-lived charged particles. 

Possible systematic biases of these limits have been considered. The relative statistical 
errors on the selection and trigger efficiencies propagate at the second order to the numer- 



ical value of the final limit p2| ; this is expected to raise by no more than a fraction of a 
percent the upper limit of the cross- sect ions attainable and, for that reason, those errors 
were not taken into account further. The efficiencies found in the search for kinks can be 
slightly overestimated for AM 11 smaller than those used in the full simulation, because 
in that case the mean momentum, and so the detectability, is lower. However, this would 
not affect the overall limits obtained, as this region is fully covered by the search for stable 
particles. Instead, at larger AM 1 * 1 the efficiencies can be slightly underestimated; in that 
case the limits displayed in figure 12 are conservative (that is the confidence level is 95% 



or higher). The same holds for the trigger efficiency. In conclusion, the limits obtained 
are reliable, or conservative, in the whole space of the SUSY parameters spanned by the 
present search. 



5 Search for charginos with ISR photons 

To look for short-lived charginos close in mass to the lightest neutralino, events with 
a few low energy particles accompanied by an ISR photon at high transverse momentum 
were searched for. Samples of e + e~ — > xfxi events with initial (and final) state radiation 
were simulated with AM ± = 0.3, 0.5, 1 and 3 GeV/c 2 . The criteria used to search for 
such events in the data were then defined on the basis of these samples and the simulated 
background samples. 



5.1 Charged and neutral particle selection 

Tracks reconstructed in the detector were considered as charged particles, to be used 
in the subsequent analysis, if their momenta were above 100 MeV/c and known with an 
error below 100%, their impact parameters were below 4 cm in the transverse plane and 
below 10 cm in the longitudinal direction and if the tracks were at least 30 cm long. 

To be accepted, neutral particles had to have an energy of at least 500 MeV. Both 
calorimetric showers and clear photon conversions in the material in front of the HPC were 



considered. Photon conversions were identified using the standard DELPHI algorithm [12 
which looks for a pair of tracks originating from a common secondary vertex and with an 
invariant mass compatible with zero. 
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Subdetector dependent criteria were used to reduce the rate of spurious neutral par- 
ticles reconstructed from electronic noise in the calorimeters or, in the case of the HPC, 
from a-emission by radionuclides embedded in the lead of the converter material. 

In the showers reconstructed by the HPC at least three of the nine layers must have 
given a signal; the first layer with a signal must be before the sixth HPC layer; not more 
than 90% of the total energy of the shower must be deposited in a single layer; the polar 
angle of the shower axis must point towards the main event vertex within 15°. 

In the FEMC a signal was required in at least two towers. In the data collected since 
1997 a more refined quality requirement was used, defining a frame of 3 x 3 glasses, 
centered on the barycentre of the shower: a shower was discarded if its energy was above 
8 GeV and more than 94% of this energy was deposited in the central glass, or if it had 
no more than three glasses hit and these were all lined up in the same row or column. 

In the STIC an energy of more than 2 GeV must be associated to the shower and at 
least two towers must be hit. 

All neutral showers inside a cone with half opening angle of 5° were combined. The 
resulting shower was not considered an ISR photon candidate if the fraction of energy 
detected in the hadron calorimeter was above 10% of the total shower energy. 

Given the reduced reconstruction and association efficiency of tracks in the two endcaps 
as compared to the barrel, spectator electrons from two-photon events can be seen as 
showers in the calorimeters without the corresponding track elements, thus faking neutral 
particles. Such electrons constitute a serious background to ISR photons, and care must 
be taken to reject these events. 

5.2 Event selection 

The analyses of the data collected before year 1997 and in 1997 have been done sepa- 
rately, with slightly different selection criteria. The tighter preselection cuts of 1997 were 
intended to reduce further the overall volume of the data used in the analysis, given the 
increased luminosity collected with respect to the previous years. The final selection of 
1997 data also exploited the improved charged particle rejection in the forward regions 
which is obtained by trying to associate hits in the Silicon Tracker to the neutral showers. 

5.2.1 Preselection 

In the preselection, events were required to have at least two and at most ten charged 
particles. There should be an isolated photon candidate of at least 4 GeV, having a 
transverse energy above 2 GeV (above 4.5 GeV in the 1997 data at the centre-of-mass 
energy of 183 GeV) and a mass recoiling against it (M opp , defined by M^ pp = E^ ms — 
2E cms E 1 ) of at least 90 GeV/c 2 . This photon had to be isolated from any other charged 
or neutral particle in the event by 15° or more. The visible energy of all particles excluding 
the photon must not exceed 8% of the available centre-of-mass energy. If all the particles 
in the event were in the same hemisphere inside a cone centered on the beam axis and with 
half opening angle of 60°, the event was discarded; this reduces the background coming 
from beam-gas or beam-wall interactions. To reject most of the two-photon background, 
also the fraction of the total visible energy within 30° of the beam axis (excluding the 
ISR photon candidate when inside that cone) was required to be less than 60%. 
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5.2.2 Selection 

Following the preselection, in order to optimize signal efficiency and background rejec- 
tion, more stringent selection criteria were imposed as follows. 

• There must be at least two and at most six accepted charged particles and, in any 
case, not more than ten tracks in the event. 

• To reject the bulk of the two-photon background, the transverse energy of the ISR 
photon was required to be greater than (K^) min , defined in eq. (jl]). Here m \ n was 
taken as the minimum angle at which the photon shower is fully reconstructed in the 
STIC (1.82°). 

• M opp must be above 96 GeV/c 2 . This was intended mainly to reduce the number of 
events with an on-shell Z recoiling against the photon. 

• The photon was required to be isolated by at least 30° with respect to any other 
charged or neutral particle in the event. 

• The sum of the energies of the particles emitted within 30° of the beam axis (£30) 
was required to be less than 25% of the total visible energy. The photon was not 
considered in any of the two energy sums if its direction was inside that cone. 

• If the ISR photon candidate was detected in the STIC, it must not be correlated 
with a signal in the Veto Counters. 

• In the data collected during 1997, if the ISR photon candidate was at an angle 
between 10° and 25° from the beam, the region where the TPC cannot be used in 
the tracking, it must not be correlated with hits in the Silicon Tracker. 

• The visible energy of the event, excluding the photon, must be below 5% of E CUIS . 
For AM ± < 1 GeV/c 2 this fraction was reduced to 2%. 

Figures §, || and [T^ show the distributions of some of the variables used for the selection 
in the preselected samples at 161, 172 and 183 GeV respectively. The data are compared 
with the SM expectation, normalized to the same luminosity. The agreement is good 
in almost all distributions. The simulated two-photon interactions giving hadrons via 
the QPM process with and without transverse momentum (jpx) of the ISR photon were 
compared. It was evident that the small excess of real data in the first bins of the 
distributions of -E30 and (E vis — Ej) can be, at least qualitatively, explained by the lack 
of high px ISR in some of the simulated two-photon samples. At present there are no 
generators available, however, which correctly describe the ISR transverse momentum for 
these processes. 

The corresponding distributions for the signal with M-+ = 50 GeV/c 2 and AM = 1 
GeV/c 2 , taken as an example, are shown to the right in the same figures. The histograms 
of the visible energy are shown for the three mass differences of 0.3, 1 and 3 GeV/c 2 , since 
the energy of the visible decay products depends on AM^. 

For the background estimates at yfs = 130 and 136 GeV, the problems with the 
simulated samples were more important: many of them were originally generated with 
insufficient statistics or with stricter requirements than those used in the subsequent 
analysis. No detailed comparison of the distributions of the data with respect to the 
simulation was meaningful at those energies. 

5.3 Results 

The results of the selection, when applied to data and simulated background at all the 
centre-of-mass energies, are shown in table [|. Six candidate events remain in the data 
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Figure 8: Some of the variables used in the selection at 161 GeV. In the left plots the 
data (dots) are compared with the SM expectations. On the right, as an example, the 
corresponding distributions (with arbitrary normalization) are shown for the signal with 
M-+ = 50 GeV/c 2 and AM ± = 1 GeV/c 2 . In the plot of the visible energy (second row) all 

three mass splittings are shown: dotted, AM ± = 0.3 GeV/c 2 ; dashed, AM ± = 1 GeV/c 2 ; 
full line, AM ± = 3 GeV/c 2 . 
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Figure 9: Same as figure but data and simulation refer to the samples at 172 GeV. 
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Figure 10: Same as figure |8|, but data and simulation refer to the samples at 183 GeV. 
The cut on the transverse energy of the photon in the preselection is tighter than the one 
used in the previous plots at lower centre-of-mass energies. 
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after the selection for 1 < AM ± < 3 GeV/c 2 (one in the sample at 136 GeV collected in 
1995, one at 172 GeV and four at 183 GeV). Some of their properties are summarized in 
table |2|. Only two of them pass also the stricter requirements for smaller AM (both at 
183 GeV). There is no significant excess above the SM expectations in either selection. 
The SM background remaining at the end of the selection is almost entirely composed of 
two-photon interactions. 



E v i s — E 
(1 < AM 

Data 


± < 3 GeV/c 2 ) 
S backgrounds 


-E v i s — £ 

(AM 1 
Data 


7 < 2% ■ -Ecms 

< 1 GeV/c 2 ) 
S backgrounds 


E cms = 130/136 GeV ( / £ = 11.7 pb" 1 ) 


1 


0.84 ± 0.84 





~ 


-Ecms = 161 GeV ( JC = 9.7 pb' 1 ) 





1.12 ±0.38 





0.45 ±0.21 


E cms = 172 GeV (J £ = 9.9 pb" 1 ) 


1 


0.64 ±0.18 





0.11 ±0.06 


E cms = 183 GeV ( J* £ = 50.0 pb" 1 ) 


4 


2.96 ±0.88 


2 


0.44 ±0.21 



Table 1: Results of the selection on the data and on the sum of the expected SM 
backgrounds. The integrated luminosity is the one used for the analysis. 

At 130/136 GeV some of the background samples were either missing or had require- 
ments at the generation level stricter than those used later in the analysis and the back- 
ground quoted in the first line of table [I] is likely to be underestimated. 

The selection efficiency (e sel ) for charged higgsinos and gauginos at yfs = 183 GeV, 



determined by using the samples of simulated events is shown in figure |TT] as a function 
of the mass of the chargino and of AM 1 . Similar efficiencies have been obtained at the 
other centre-of-mass energies. It must be stressed that the very low signal efficiency comes 
from the requirement of having an energetic ISR photon radiated at visible angles. As 
the mass of the chargino increases, the energy of the photon decreases, thus lowering 
the overall event selection efficiency. The selection efficiency for gauginos, in the case 
of a heavy sneutrino, is higher than for higgsinos, even though the decays are similar 
in the two scenarios, because the gaugino cross-section resonates more strongly around 
the Z pole, leading to more frequent ISR radiation. However, the gaugino efficiency is 
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2 


1 
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2 


1 
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2 


2 
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0.05 ± 0.04 
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183 


2 


1 


13.2 


5.8 
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0.00 ±0.01 


0.00 ±0.01 



Table 2: Some of the properties of the events remaining in the data after the selection: 
the centre-of-mass energy, the number of charged and neutral particles in the event, the 
energy and the polar angle of the photon, the impact parameters in the plane r0 of the 
two charged tracks. 



significantly smaller than the one shown in the figure when the mass of the sneutrino is 
close to the mass of the chargino: light sneutrinos would enhance the fraction of leptonic 
decays, which have additional missing energy and then a lower efficiency for the same 
AM 1 * 1 . All this was taken into account when computing the exclusion limits. 

Below transverse momenta of 1-2 GeV/c the single track trigger efficiency of DELPHI 
starts to decrease PU[ . The events searched for have few tracks of low energy (besides 
the energetic photon); for that reason, a dedicated and detailed study of the trigger 
performances at low visible transverse energy was carried out. The trigger efficiencies for 
events with a visible transverse energy similar to the one of the system accompanying the 
ISR photon in the signal samples which pass the selection cuts, were computed from all 
the data collected at LEP2. This study had to be repeated for every new year of data 
taking, since the definitions of the decision functions and the downscaling factors applied 
in case of noise during the data acquisition can affect the results. The sample of events 
with independent trigger signals from any of the calorimeters was used to compute the 
efficiency for the tracking components of the trigger. The overall trigger efficiency for the 
signal (e tr& ) was then estimated as the logical OR of the single photon trigger efficiency 
23|j and the trigger efficiency for events in which there are only low energy tracks. For 



the signal events considered here, the estimated overall trigger efficiency varies between 
82% and 98%. 



5.4 Exclusion limits in the search with the ISR tag 

Since there was no evidence for a signal at any of the centre-of-mass energies studied, 
limits were derived on the pair production of charginos nearly mass-degenerate with the 
LSP. 

5.4.1 Method 

Although the cross-section is relatively high for all chargino masses (see figure ^J), it is 
important to combine the data taken at all energies, because of the low signal efficiency. 
The method used to obtain a combined limit from the data taken at the different centre- 
of-mass energies is the following. 
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Figure 11: Selection efficiencies for charged higgsinos and gauginos at 183 GeV as a 
function of their mass and mass difference with the lightest neutralino. This gaugino 
efficiency is valid only in the approximation of heavy sneutrinos; for light sneutrinos a 
different value is used. 



The number of candidates observed at each centre-of-mass energy is distributed ac- 
cording to a Poissonian. The (Bayesian) probability density of the true mean value for 



the number of signal events in the data (iV si g) is derived according to [24], given the num- 
ber of observed events and the expected number of background events at each -E C ms- The 
error on the expected background content was taken into account by assigning a Bayesian 
probability density to different values of the background and weighting the iV S ig densities 
by the background density. The background densities used were taken to be Gaussian at 
every E cmB , with standard deviations equal to the errors reported in table |l|. 
The total number of expected signal events is 

N eKp = 1 £<T i £ i e"*e it * (4) 

where i runs over all the centre-of-mass energies. The statistical errors on the selection 
efficiencies completely dominate over the errors on the integrated luminosities and trigger 
efficiencies, which were neglected. For the cross-sections, the lowest values obtained in the 
scan of the SUSY parameters were used. iV ex p was then assigned a Bayesian probability 
density, based on the binomial statistics relevant in the calculation of e sel , and assuming 
equal a priori probabilities. 

The two probability densities for iV S ig and iV ex p are independently determined, and the 
probability that iV S ig < iV ex p can therefore be obtained by convoluting these two densities 
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using Monte-Carlo techniques. If this probability, in a given scenario, is equal to r\ for a 
point in the plane (M$+, AM ± ), then the point is excluded at the confidence level r\. 

In this way, confidence levels of exclusion were derived for the mass points where a 
full simulation of the signal had been performed. Between these points an interpolation, 
based on SUSYGEN events without full detector simulation, was used to obtain the limit. 

5.4.2 Limits 

Table § gives the 95% CL lower limits on the mass of the chargino for each of the three 
scenarios considered and for the different AM^ ranges. To compute the efficiencies, event 
samples have been generated at fixed values of the mass of the chargino and of AM . 
While it is straightforward to interpolate between simulated points at the same AM ± and 
different M~+, it is more difficult in the case of different DMP because of the different 
Q-values in the chargino decay. For this reason, a limit was calculated for each AM^ 

simulated by interpolation in M-+, but no interpolation in AM^ was done. Instead, 

Xi 

the limit between any two simulated AM^ was conservatively taken as the the lower of 
the corresponding two M-+ limits, giving a step-like exclusion contour. In particular, 
no limit was derived below the minimum AM^ used in the simulation; however, using a 
small sample of simulated events it has been verified that those efficiencies drop quickly 
below that lowest AM ± , at the centre-of-mass energies studied in this paper. 
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62.6 GeV/c 2 
62.6 GeV/c 2 
60.6 GeV/c 2 


49.4 GeV/c 2 
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Table 3: 95% CL lower limits on the mass of the chargino obtained with the search for 
soft particles accompanied by an energetic ISR photon, in the three scenarios in which a 
mass-degeneracy with the neutralino is possible. 

When the chargino is a gaugino and the mass of the sneutrino is above but close to 
M-+ , the decay xt ~^ X^ +U i > mediated by a virtual v, is enhanced with respect to the 
branching ratios shown in figure |3|. In that case there are two invisible particles which 
carry away the energy and the visible lepton is usually softer than the pion of the two- 
body decay xf ~* Xi 1T+ - F° r this reason the selection efficiency becomes lower with a 
light sneutrino with mass above the chargino mass and the limit has been set only down 
to AM ± = 500 MeV/c 2 in table |3|. If the sneutrino is lighter than the chargino, whether 
or not it is the LSP, the AM to be considered is the mass difference between the chargino 
and the sneutrino: if M„ > M^o the subsequent decay v — > vx\ is experimentally invisible; 
if < M^o the chargino nevertheless decays predominantly into Iv. 

All these results take into account a variation of tan f3 between 1 and 50 and a variation 
of the Mi, M 2 and /x parameters so that the mass difference between the chargino and 
the neutralino remains below 3 GeV/c 2 and M 2 < 2M X < 10M 2 (although a coarse scan 
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of the space of the SUSY parameters indicates that all limits would remain the same even 
HM 1 >5M 2 ). 

As anticipated in section |2|, the limitations of some of the two-photon generators used 
lead to an underestimate of the SM background in the region of high ISR where 
the signal is expected. There is no evidence of any signal, however, and when deriving 
exclusion limits the only bias expected because of the somewhat inadequate simulation 
is that the mass limits are likely to be underestimated. In other words, the limits are 
conservative, in the sense that with a more precise background simulation their confidence 
level would be likely to exceed 95%. No attempt to compensate for this effect has been 
made. 

Figure |12| shows these limits together with the ones obtained in the search for long- 
lived charginos. For completeness, also shown are the limits obtained at LEP1 |5||| and 
the results of the search for high AM ± charginos in DELPHI B . 



6 Conclusions 

Charginos nearly mass-degenerate with the lightest neutralino (assumed to be the LSP) 
have been searched for using the data collected by the DELPHI experiment during the 
runs of LEP above the Z pole energy from 1995 to 1997. Two different approaches were 
used. 

For extremely small mass differences (AM < 200 MeV/c 2 ) the information contained 
in the long lifetime of the chargino (decay length) was exploited. The search for heavy 
particles decaying outside the central subdetectors of DELPHI used the specific ionization 
of the tracks in the TPC and the light cone produced in the RICH detector. Chargino 
decays, between about 15 cm and 1 m, would be seen as kinks in the tracks reconstructed 
in the central tracking detectors of DELPHI. No candidate events were found in the data 
collected at 130, 136, 161, 172 and 183 GeV using the first analysis, or at 161, 172 and 
183 GeV using the second analysis. 

When the lifetime is so short that the decay vertex cannot be seen inside the sensitive 

devices of DELPHI and AM ± is too small (0.3 < AM* < 3 GeV/c 2 ) to be selected by 

the usual criteria adopted in the search for charginos, some events can still be recovered 

by looking for the typical topologies of the chargino decays at low AM* accompanied 

by a high energy photon radiated from the initial state. The ISR signature reduces 

the otherwise overwhelming two-photon background to acceptable rates, although it also 

strongly affects the signal efficiency It is necessary to combine all the statistics collected 

so far at the different LEP2 centre-of-mass energies to achieve sufficient sensitivity. No 

evidence of a signal has been found in the data collected by DELPHI at the centre-of-mass 

energies of 130, 136, 161, 172 and 183 GeV. 

The regions of the plane M-+ vs. AM ± excluded at the 95% CL by the combination 

*i 

of these searches in DELPHI are summarized in figure [12]. There is still an inaccesible 
region at, approximately, 200 MeV/c 2 < AM ± < 300 MeV/c 2 . With the higher statistics 
available and the increased boost of the decay products of the chargino, provided by the 
raised centre-of-mass energy, this region is likely to be explored with the use of data from 
1998 onwards. 
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Figure 12: Regions in the plane (M-+,AM ± ) excluded by DELPHI at the 95% CL 

using: the search for high AM 1 * 1 charginos; the search for soft particles accompanied by 
ISR; the search for long-lived charginos. The three scenarios are the ones which allow low 
/S.M ± : the lightest chargino is a higgsino; the lightest chargino is a gaugino and Mp > 500 
GeV/c 2 ; the lightest chargino is a gaugino and M v > M-+ (or M v < M- 
replacing AM ± on the y axis in the search with the ISR signature). 
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